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SUMMARY 

The constant magnetostriction coefficient for nonferromagnetic materials is derived 
in terms of changes in the Fermi surface. The theory is applied to Fermi surfaces of 
spherical and elliptical geometries. The magnetostriction coefficient along the trigonal 
axis of bismuth is calculated and compared with experimental results. 


INTRODUCTION 

The magnetostriction of nonferromagnetic materials was first observed by Kapitza 
(ref. 1) for large transient fields. Subsequently, Shoenberg (ref. 2), Wolf and Goetz 
(ref. 3), Anderholm (ref. 4), and Aron (unpublished NASA data) investigated this phenom- 
enon. These investigators observed that magnetostriction was a function of crystal direc- 
tion relative to the field and was dependent on the first power of the magnetic field. The 
phenomenon was assumed to arise from the change of magnetization, which provides the 
energy for elastic deformation. 

The magnetostriction coefficient of a sample in a homogeneous magnetic field is de- 
fined as the fractional change in length per unit magnetic field. The magnetostriction is 
called longitudinal or transverse depending on whether the change in length is measured 
in the direction of the magnetic field or perpendicular to it. 

In classical physics, a system of free electrons confined to a fixed volume has zero 
diamagnetic moment. In quantum mechanics, the quantization of angular momentum in a 
magnetic field results in a nonzero magnetic moment. The magnetization is composed of 
two components, one of which is constant and the other is oscillatory with respect to 
changes in magnetic field. This report is concerned only with the constant component of 
the diamagnetic moment. 

The constant component of magnetization is formed chiefly from the bunching of en- 
ergy levels. The application of a magnetic field does not, however, change the average 
density of the levels. The magnetostriction coefficient, which depends on the magnetiza- 
tion, is derived in terms of the effect of an elastic deformation on the extremal area of 



the Fermi surface. This deformation is caused by a change in both the equilibrium elec- 
tron density and the energy of the electrons. 

hi general, an electron in a metal is governed by a complicated dispersion law. The 
frequently used quadratic dispersion law for conduction electrons is valid only when the 
Fermi surface is the lower part of the Brillouin zone. The theory presented in this re- 
port can be applied to an arbitrary dispersion law. 
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longitudinal magnetostriction 
speed of light 
electron energy 

electron energy in undeformed state 
electron energy in deformed state 
electron charge 

Fermi distribution function 

density of states 
magnetic field 

Plank's constant divided by 2 ir 
compressibility (isotropic), 
deformation parameter 
Boltzmann constant 
wave number 

normal displacement of Fermi surface in k space 

radius of spherical Fermi surface 

deformed length 

undeformed length 

magnetization per unit volume 

temperature component of magnetization 

effective mass tensor element 
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mass of free electron 

number of electrons per unit volume 

number of electrons per unit volume in undeformed state 
pressure 

electron momentum 
extremal momentum 
electron momentum in z-direction 
summation index 

cross-sectional Fermi surface area 

compliance tension element 

extremal cross-sectional area of Fermi surface 

temperature 

volume of sample 

Fermi velocity 

m 0 / m ij 

hole effective mass 
phase 

directional cosines 
dilation 

see equation (8) 

Kronecker delta function 
chemical potential 

chemical potential in undeformed state 
stress tensor element 
strain tensor element 


THEORETICAL ANALYSIS 

The longitudinal magnetostriction coefficient developed from thermodynamics (ref. 5) 
is given by 
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The first term of equation (1) gives an extremely small contribution (about 10” ) per 

gauss. Therefore, the magnetostriction coefficient is very closely given by 


C = - 



( 2 ) 


The constant diamagnetic component of the magnetization of an electron gas as given 
by Lifshitz and Kosevich (ref. 6) is 



where 

f (x) s (l + e x ) 1 

and 3S/3p 0 is a derivative of S with respect to p evaluated at the extremal momentum. 
For kT « p, equation (3) becomes 



For T » 0, the second term in brackets contributes very little, and therefore, 
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If a function cp is considered, where 
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and equation (4b) becomes, 


C = - 



R=°° 
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To determine the component dp/da, the following reasoning is applied. The magnetic 
field disturbs the lattice, moving the atoms away from the position of zero field equilib- 
rium and thereby deforming the Fermi surface and changing the value of p. It is assumed 
that the influence of an elastic deformation of the latter can be taken into account in the 
form of a small addition to the electron energy proportional to the deformation tensor. 

This assumption implies that the energy of the electrons after deformation is 


*. 


E(k) = E 0 (k) + E’(k) (6) 

where E(K), E q (K) are the electron energies in the deformed and undeformed metal, re- 
spectively, and E’(k) is the change in the energy of the electrons due to a deformation. 

Assuming the number of electrons is constant and using equation (6), Ziman (ref. 7) 
showed that the change in the chemical potential is 
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The changes involved in the energy surfaces may be described by a deformation ten- 
sor K^(k) introduced by Pippard (ref. 8). This tensor is defined by 
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( 8 ) 


“n="iiV k > 

where is the strain tensor element and 6k n is the normal displacement of the 
Fermi surface. The right side of equation (8) represents the sum of the nine terms that 
results from summing the dummy indices i and j over the range of values 1, 2, and 3. 
Under a deformation, the energy of an electron is changed by the amount 

E ’ = [ v k ( E) ] 5k n (9 > 

With equation (8), E’ can be written as 

E T = hVp^K.^k)^.. (10) 

Substitution of equation (10) into equation (7) leads to 
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If the generalized form of Hook’s law = S ijk£ a k£ is used, where is the 

elastic compliance tensor, the change in chemical potential with stress can be written as 
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Hence, from equation (5), 
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For temperatures above T = 0, the second term in equation (4a) must be included, 
and this gives the contribution 
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Hence, the temperature component of magnetostriction coefficient can be written as 
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APPLICATIONS 

The longitudinal magnetostriction coefficient is computed for two cases. As an ex- 
ample of a spherical Fermi surface, sodium is used, and for other alkali metals, y= 1/2; 
hence, 
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For a spherical Fermi surface, 


K.. = - 6..k 
i] 3 H o 


where k Q is the radius of the Fermi surface. 

Inserting the previous values into equation (13) yields 
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( 21 ) 


where 


k = 156 cirr /dyne 


p = 3.2 eV 
o 


therefore, 


C/H = 1.85X10 -18 , 


which is too small to be observed experimentally. 


ELLIPSOIDAL FERMI SURFACE FOR BISMUTH 

Before calculating the coefficient of magnetostriction, it is advantageous to explain 
some of the salient features of bismuth. Bismuth has five electrons per atom and is 
classified as a semimetal. This metal has a small number of free electrons due to over- 
lap outside the Brillouin zone. It is believed that the Fermi surface of bismuth consists 
of highly anisotropic ellipsoids. 

No detailed calculation of the energy bands of bismuth has been made. Thus, it is 
impossible to determine K„(k) and effective mass analytically, but it may be possible 
to evaluate them by experimental techniques. Fortunately, a large amount of experimen- 
tal information has accumulated on the band structure of bismuth. Measurements of the 
de Haas - Van Alphen cyclotron resonance and ultrasonic attenuation can be interpreted 
in terms of a simple model of ellipsoidal energy surfaces. An interpretation of the 
de Haas - Van Alphen experiments for electrons has shown that the axis of an ellipsoid is 
not along but is tipped out of the principle symmetry axis of the crystal. Other models 
propose more than three electron ellipsoids; usually they assume six ellipsoids. Most of 
these models have one factor in common, which is that some of the ellipsoids are asso- 
ciated with holes and some with electrons. The question of the correct model for bismuth 
has not yet been resolved. 
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Recent investigations by Aron (unpublished data), and Brandt and Venttsel (ref. 9) 
have indicated that some of the physical properties for bismuth, such as the amplitudes 
on the trigonal direction of oscillatory magnetostriction and very low temperature 
de Haas - Van Alphen effects, come from the hole ellipsoid. 

By measuring electrical resistivity and the Hall coefficient, Jain and Jagge (ref. 10) 
determined that hV(k)K^(k) along the trigonal axis is about 2. 5 electronvolts. 

In order to calculate the magnetostriction coefficient, two models are employed. For 
the first model the following assumptions are made: 

(1) The electron ellipsoids contribute negligibly to magnetostriction. 

(2) The term KV(k)K^(k) is due only to the hole contribution. 

(3) The hole contribution consists of two identical hole ellipsoids. 

When the constant energy surfaces are no longer spherical, a scalar effective mass does 
not apply. Since equation (19) can no longer be used, equation (13) must be referred to. 
The hole ellipsoid is described by the dispersion law 

2m o?h “ % P x + %2 p ? + Vz < 22 > 


The extremal area S m (pj i ) of cross section of the ellipsoid can be evaluated by ana- 
lytic geometry from equation (22) to be 


S (p) = 
m VK ' 
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where 


X h = r ll^ /3 22 /3 33^ + r 22^ /3 ll /3 33^ + r 33^11^22^ 

and are the directional cosines that depend on the orientation of the magnetic field, 

and the subscripts on X refer to the holes. In the following, the theory of the case where 
the field is parallel to the trigonal axis is applied. 

The area of the curve of constant energy is obtained by using the Bohr - Sommerfeld 
quantization condition and can be shown to be 

S(p„,P z > =//d P x dP y =(n + i)M£H= S m (p h ) 
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Hence, the following can be obtained: 
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From quantum mechanical considerations y is restricted to the range 0 <y < 1. 
For the case of a quadratic dispersion law y = 1/2, and for other Fermi surfaces y dif- 
fers from 1/2. 

Using the experimental data obtained by Jaggi and Jain (ref. 10) for the term 
hV(k)K^(k) and substituting equation (24) into equation (13) as well as summing over the 
two ellipsoids yield 
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For the second calculation, the nontilted eight ellipsoid model proposed by Abeles 
and Meiboom (ref. 11) will be used. This model proposes two hole ellipsoids and six 
electron ellipsoids. For this case, the assumption employed will be that the holes and 
electrons contribute equally to the term KV(k)K..(k). Using equations (21) to (24) and 
substituting p g and produces the electron's Fermi energy and the electron's effec- 
tive mass for p h and /3y, respectively. Summing over the eight electron ellipsoids 
yields the following: 
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To obtain the two hole ellipsoid contribution equation (25) is employed. The magneto- 
striction coefficient for this model is, therefore, 
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A factor of 2 appearing in the term [hVp(k)K-(S)S^j^]/2 is due to the assumption that 
the electrons and holes contribute equally. 

From cyclotron resonance, de Haas - Van Alphen, and ultrasonic attenuation data 
(ref. 12), the following values are obtained: 


On = 2. 02 

a 22 = 67 

a 33 = 83 ' 3 
P e = 2. 9X10” erg 

Ph = 2.6xl0“14 erg 

S ijk£ = 4x10 12 cm 2 /dyne 

K = 3. 2X10" 12 cm 2 /dyne 
J8 n = 2.0 

/ 3 22 = 2.0 
/3 33 = 1.4 

= 1. 76X10" 16 

/trig 

— (electron + holes) = 1. 86x10" l 8 
H 

The experimental values of C/H are reported in table I. 

TABLE I. - RATIO OF LONGITUDINAL 
MAGNETOSTRICTION TO MAGNETIC 


FIELD IN TRIGONAL DIRECTION 


Source 

Trigonal values 


of C/H 

Wolf and Goetz (ref. 3) 

5.7XHT 16 

Kapitza (ref. 4) 

6.5 

Anderholm (ref. 4) 

8.3 

Aron (unpublished NASA data) 

7.2 
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CONCLUSIONS 


The theory presented for the magnetostriction coefficient cannot be accurately tested 
for bismuth. In order for the theory to be adequately employed, a band calculation for 
this metal must be made. An alternate procedure would be to obtain sufficient experimen 
tal data on the deformation parameter for bismuth. 

If this theory is correct, the value of the magnetostriction coefficient C depends 
basically on two quantities: the change of Fermi surface area with momentum, and the 
deformation parameter. If the deformation parameter K„(k) can be evaluated from 
theoretical calculations and C can be determined from experimental measurements, the 
phase factor y maybe determined for complicated surfaces. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, May 12, 1966, 

129-02-05-09-22. 
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